Ground State Properties of Cold Bosonic Atoms At Large Scattering Lengths 
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In this Letter, we study bosonic atoms at large scattering lengths using a variational method where 
the condensate amplitude is a variational parameter. We further examine momentum distribution 
functions, chemical potentials, the speed of sound, and spatial density profiles of cold bosonic atoms 
in a trap in this limit. The later two properties turn out to bear similarities of those of Fermi gases. 
The estimates obtained here are applicable near Feshbach resonances, particularly when the fraction 
of atoms forming three-body structures is small and can be tested in future cold atom experiments. 



Bosc-Einstcin condensates(BECs) near Feshbach res- 
onances have been one of the most exciting ultracold 

f stems studied so far in experiments [l|, 0, y|, EL 0, S 
H, @| . On one side of the resonances where the scat- 
tering lengths are negative, fascinating collapse-growth 
cycles due to thermal clouds and spectacular controlled 
collapsing-exploding dynamics have been observed [H 0, 
3] and studied theoretically On the other side, to- 
wards the resonances where scattering lengths are pos- 
itive, strongly repulsive ultracold bosonic atoms and 
their intriguing properties have been explored 0, @, 0, HI]. 
Despite the reduced lifetime of the cold gases in this 
limit due to enhanced three-body recombination, quite 
remarkable progress has been made to probe interac- 
tions between atoms. Recently, pursuit in this direc- 
tion has been revived, and more vigorous efforts have 
been madefl, Our theoretical studies in this Letter 
are mainly motivated by these experiments. Bose gases 
at large positive scattering lengths have been a horren- 
dously challenging topic in theoretical physics for more 
than half a century [3, EH, El E! E3- The standard 
low density expansion that works quite well for dilute 
gases is not applicable when the scattering length a is 
comparable to or even larger than the mean interparti- 
cle distance d. Here we suggest a variational approach 
which takes into account two-body correlations and can 
be extended to the limit of a large positive scattering 
length. We further apply this approach to estimate var- 
ious fundamental properties of cold bosonic atoms near 
Feshbach resonances, particularly when the fraction of 
atoms forming three-body structures is small. Unique 
features in the momentum distribution function, chemi- 
cal potential, speed of sound and the cold atom density 
profile in a trap can be potentially probed in experiments. 

Cold bosonic atoms at large scattering lengths were 
also previously addressed in a few inspiring theoretical 
papers [H, Eil- Cowell et al. estimated chemical po- 
tentials and condensate fraction s by employing distinctly 
different Jastrow wavefunctions[15l| . There are a few in- 
terconnected differences between their results and ours. 
First, while the physics at distances much shorter than 
the mean interparticle distance d is described quite ac- 
curately by the Jastrow wave functions, basic aspects of 
the long wavelength physics are not expected to be well 



captured. On the other hand, our wave function is con- 
structed under a constraint in Eq. and captures essen- 
tial features of low energy collective properties of BECs. 
For instance, the momentum distribution function has 
a -r divergence near k = for all scattering lengths, and 
at short distances our wave function is almost identical to 
the solution to the Schrodinger equation for two interact- 
ing atoms. Second, since the contribution to the deple- 
tion fraction, or the fraction of atoms occupying nonzero 
momentum states, is mainly from low energy states, we 
expect that our results are more reliable. In fact, we find 
that the depletion fraction reaches a constant value of 
about 0.5 near resonances. On the contrary, the conden- 
sate fraction estimated in Ref.pjl quickly reaches zero 
when the scattering length a becomes comparable to d, 
suggesting that atoms could be completely depleted from 
the zero momentum state and there should be an unex- 
pected quantum phase transition at a finite scattering 
length. Third, chemical potentials near resonances esti- 
mated there appear to be bigger than the values obtained 
in our calculations. This seems to imply that the trial 
wave functions adopted here should be an energetically 
better candidate for ground states. 

The trial wave function in Eq.® effectively encodes 
two-body correlations. To include high-order correla- 
tions such as three-body effects, a much more sophistic 
ansatz is needed. A nontrivial role of three-body inter- 
actions was previously appreciated by Braaten et al. in 
Ref. where the effects on BECs were estimated in 
the limit of a small scattering length. Although there 
was no definite evidence of Efimov trimers in BECs of 
sodium or rubidium atoms studied in Ref. [H, @, 0, 0] , an 
earlier experiment on relaxation rates of cesium atoms 
did show, as a precursor of two-body resonances, addi- 
tional structures which had been attributed to Efimov 
states [13, El, El [H. More efforts are to be made to 



understand the nature of BECs in this limit and the ap- 
proach proposed below is a baby step towards this di- 
rection. Our results are valid when the three-body cor- 
relations induced by Efimov trimers are not dominating. 
The question of whether the emergence of Efimov trimers 
introduces distinct modulations to the scaling functions 
discussed below, or mainly sets the lifetime of BECs, rep- 
resents an exciting new direction that is worth pursuing. 
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Moreover, our scaling hypothesis works best when the 
typical range of interactions tq is much less than the 
inter-particle distance d. When the density increases, de- 
viations from the scaling behaviors become substantial, 
and the scaling functions proposed below are no longer 
sufficient for characterizing BECs. Eventually, a quan- 
tum gas might undergo a transition to a dense liquid 
phase when ro becomes comparable to d. For cold atoms, 
this fortunately only occurs at a density which is not ex- 
perimentally accessible because of severe trap losses. 

We consider bosonic atoms that interact with a short 
range potential of range r and scatter at two-body scat- 
tering lengths <z(> 0). For BECs with a number density 
poi assuming two-body effects are dominating, we can 
generally express the momentum distribution function 
n k and the chemical potential /i in terms of dimensionless 
functions / and h, i.e., n k = f(kd, §, ^), fi = e F /i(§, ^); 

and d = (i- 3 — ) 1 ^ 3 , £f = ^ 67r „ Po ^ — . For short range inter- 
actions, ro is much smaller than the mean inter-particle 
distance d so that we approximate ^ to be zero but a 
can vary over a range from much smaller than d to much 
bigger than d. Function / and h thus depend only on 
two dimensionless variables x — kd and y — § and are 
reduced to two scaling functions f(x, y) and h{y), respec- 
tively. The functional form of f(x,y) and h(y) proposed 
in this way does not depend on details of interaction po- 
tentials or number densities or scattering lengths and is 
universal; / and h characterize basic properties of BECs. 
Note that, as illustrated below, the Fermi energy cf that 
is normally defined for a Fermi gas with the same number 
density po turns out to be the only relevant energy scale 
for BECs near resonances. 

When a is much smaller than d, these functions can be 
obtained by using the standard mean 

field [Him GSEi]. 

Indeed, in the dilute limit when y is much less than unity 
one can verify that 
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where we also introduce g(y) for the depletion fraction. 
f(x, y) is divergent as yfyjx when x or momentum k ap- 
proaches zero; this behavior is an indication of gapless 
soundlike collective excitations in BECs. Furthermore, 
that f(x,y) decays as y 2 /x 4 in the large- a; or large-fc 
limit reflects the free particle nature of high energy ex- 
citations. For cold atoms at large scattering lengths, y 
is substantial and the form of / and h functions remains 
to be understood. In the following we are going to inves- 
tigate these scaling functions in the limit when a(or y) 
becomes comparable to or bigger than <i(or 1). 

To quantitatively study / and h functions in the limit 
of a large scattering length, we adopt a variational ap- 
proach to BECs. In this method, Co, the condensate am- 



plitude and gk, k ^ 0, the pairing amplitude that is re- 
lated to the occupation number of atoms in a state of 
momentum k, are variational parameters. We then min- 
imize the energy with respect to <7 k and cq but with the 
total number of atoms Nt fixed. 

To introduce trial wave functions which are viable in 
both small and large scattering- length limits, we require 
that at any given scattering length the ground state 
should be a vacuum of Bogoliubov quasi-particles and 
is annihilated by a set of quasiparticle operators 
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|g.s.) = 0. (2) 



Here cik (a^) is an annihilation (creation) operator for 
an atom with momentum k. Detailed structures of the 
quasi-particle operators are specified by real variables <?k 
and will be determined variationally below. The ansatz 
that satisfies Eq.((2]) can be written as 



exp 



k z>0 



(co4) n exp ( 5k a k «-k) l°>- ( 3 ) 



Here A is the normalization factor. Again Co is the 
condensation amplitude and is the pairing amplitude 
with |(?k| < 1; for ground states, we further assume 
g_k = gin- This trial wave function encodes two-body 
correlations but not three-body ones. Similar wave func- 
tions have been used to study pair condensates of at- 
tractive bosons 21]. nk, the occupation number of atoms 
with momentum k, is a simple function of g^ 



n k = (a k a k ) = 



1.9k 



The Hamiltonian of cold bosons is 
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H = J2 
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J7(q)a kl a k2 , (5) 



U(<l) — hi d 3 rU(r) exp(iq • r) is a two-body interaction 
potential, and is the volume of the system. The total 
energy Et of the trial state \g.s. > is evaluated to be 
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FIG. 1: (Color online) The momentum distribution function 
n(k) for different interparticle distance d and scattering length 
a [in units of ro , the range of interaction and in a, k is in units 
of h/ro ]■ In all cases, n(fc) has a desired 1/k divergence when 
k approaches zero. In b), n(fc) functions for different d but 
with the same value of y = ^ are further shown to collapse to 
a single scaling function when plotted against x = kd. The 
resultant three curves are for f(x,y) with y — 0.1, 1 and 10 
(from bottom to top). Depletion fraction g(y) is plotted in 
the inset; the mean field(MF) g(y) in Eq.flJ is also shown as 
a reference. 



To facilitate discussions on large scattering lengths, we 
assume that the interaction potential is a square well 
one, U(r) = —V when r < ro but otherwise is zero. 
The corresponding s-wave scattering length is a — ro — 

tan (^JmUroj / \fmU . We choose the depth of the po- 
tential U to be 7r/2 < V mUro < n so that ro < a < oo. 

To obtain ground states, we minimize the total energy 
in Eq.® with respect to parameters g^ and Co, subject 
to a constraint that the total number Nt is fixed, 



N T = |c | 
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When the potential is weakly repulsive, we verify that 
the minimization does lead to the standard results for 
weakly interacting BECs, i.e. Eq.{T]). For attractive 
potentials introduced above, the minimization is carried 
out numerically. When scattering lengths are positive, 
one of the energy minima turns out to be a collection of 
molecules as expected from a two-body consideration; in 
these molecular states, the condensate amplitude is found 




FIG. 2: (Color online) Chemical potentials fi and sound ve- 
locities v s as a function of scattering length, a) is for /i (in 
units of l/2md 2 , d — 50 versus a (d and a are in units of ro); 
in b) we further plot fi in units of the Fermi energy ef as a 
function of y = a/d and illustrate two plots in a) collapse into 
a single scaling curve. The resultant plot defines the scaling 
function h(y)(= h/sf) for an arbitrary y. The dashed line is 
for h(y) in Eq.([TJ. Shown in the inset is v s (in units of the 
Fermi velocity vf ) versus a/d. 



to be zero and \g^\ is less than unity for all k. To under- 
stand BECs of scattering atoms in open or nonmolccular 
channels that are most relevant to experiments on cold 
atoms, we project away the molecular states and mini- 
mize the energy in the subspace of scattering channels. 
This is achieved by imposing a projection constraint on 
fk, X)k5k l °'5k = 0> where g^° l are the calculated values 
of gk for molecular states. This vanishing inner product 
between molecular states and states of atoms effectively 
projects out a desired subspace of open channels. 

Below we present results for BECs with different densi- 
ties and scattering lengths. The minimization algorithm 
does converge leading to a ground state in the subspace 
when we set <7k + 1 to be proportional to k in the close 



vicinity of k = 0[22|. We further find that decays as p- 
in the large-fc limit for all scattering lengths. Following 
the relation between <?k and in Eq.([4|), one then ob- 
tains the asymptotics of nk in both large-/c and small-fc 
limits. The characteristics in these two limits are ro- 
bust and, when the scattering length a is tuned, remain 
the same as those in Eq.([T]). However, the crossover en- 
ergy between these two limits, which is approximately 
the chemical potential, strongly depends on the scatter- 
ing lengths or densities (see FigQ}. When plotted against 
x = kd, data for or n(k) calculated for different den- 
sities and scattering lengths all collapse to a single set of 
curves which correspond to f(x,y) for different y = §. 
Furthermore, we observe that the function n(k) — f(x, y) 
quickly approaches foo(x) when y exceeds unity. Using 
the momentum distribution function, we also estimate 
g(y), the fraction of atoms that are depleted from the 



4 




(aN 1 7L HO ) 1 ' 3 ho 

FIG. 3: (Color online) a) Radius R of BECs in a spher- 
ical harmonic trap with harmonic length Lho as a func- 
tion of scattering length a; the dashed line is the mean 
field Thomas-Fermi radius Rtf (see discussions before Eq.© 
). b) Spatial density profiles in a harmonic trap at differ- 
ent scattering lengths; the density at the center p(0) is esti- 
mated to be p(0)a 3 = 0.00012, 0.038, 0.24, 180 for the dashed- 
dotted, dashed, dotted, and solid line respectively. In the 
inset, we also plot the mean field result(black dashed line) for 
p(0)a 3 = 180 for a comparison. In a),b), the total number of 
atoms is set to be Nt — 2 x 10 6 . 

zero momentum state; g(y) saturates at a constant value 
near resonances. 

The chemical potential is studied by evaluating /i = 
OEt/ONt- In the limit of large scattering length, the 
main characteristic is that ii saturates at a value of 
around 80% of the Fermi energy tp of the correspond- 
ing density. When the chemical potential in units of tp 
is plotted against scattering lengths y = §, all data again 
collapse to a single master curve which quantitatively de- 
fines the scaling function h(y)(= introduced above; 
and h(y) approaches 0.8 once y becomes much bigger 
than unity (see Fig[2]). v s , the speed of sound that de- 
pends on the compressibility of BECs, can also be ob- 
tained by using the general relation v 2 s = po/m(d[j,/dpo). 

The scattering-length dependence of the chemical po- 
tential discussed here implies a very peculiar evolution 
of sizes of BECs in a trap (with a harmonic length 
Lho) when scattering lengths a are increased. In the 
limit of a small scattering length the size of conden- 
sates increases as a function of scattering length a and 
the mean field Thomas-Fermi radius in a spherical trap 
is R T f/L H o = (l5N T a/L HO ) 1/5 \M- As the chemical 
potential saturates at a value of 0.8eF when scatter- 
ing lengths become much larger than the typical inter- 
particle distance in a trap, the radius of the BECs in this 
strongly interacting regime is also expected to approach 
a value of 

= AN*' 6 ; (8) 

Lho 



Numerical calculations further show that A = 1.9. As 
another application of our variational approach, we quan- 
titatively investigate radii of BECs near resonances using 
a local density approximation (see Fig|3|). 

In conclusion, we have examined basic properties of 
cold bosonic atoms at large scattering lengths. Using 
the variational method, we estimate various properties 
that can be potentially tested in future cold atom ex- 
periments. Near resonances, we have found that the 
chemical potential, speed of sound, and the spatial den- 
sity profile of cold bosons in a trap resemble the cor- 
responding properties of Fermi gases. This particular 
aspect is also a unique feature of one-dimensional Tonks- 
Girarde au g ases where bosons are viewed as fermionized 
particles [3 [H, [26}. Our results are applicable near Fes- 
hbach resonances but before the Efimov physics fully sets 
in. This work is supported by NSERC, the Canada and 
Canadian Institute for Advanced Research. We thank 
Jason T. L. Ho for a stimulating discussion. 
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